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Abstract 

We examine non-thermal dark matter production from a late-decaying scalar field, with a particular 
attention on non-renormalizable operators of Z) = 5 through which the scalar field decays into the 
standard model particles and their superpartners. We show that almost the same number of 
superparticles as that of particles are generally produced from the decay. To avoid the gravitino 
overproduction problem, the decay is favored to proceed via interactions with an intermediate 
cut-off scale M <^ Mp. This should be contrasted to the conventional scenario using the modulus 
decay. The bosonic supersymmetry partner of the axion, i.e., saxion, is proposed as a natural 
candidate for such late-decaying scalar fields. We find that a right amount of the wino/higgsino 
dark matter with a mass of O(IOO) GeV is obtained for the saxion mass around the weak scale and 
axion decay constant, F4 = O(10^~^^) GeV. 
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I. INTRODUCTION 



The lightest superparticle (LSP) is one of the most plausible candidates for the dark 
matter (DM) in supersymmetric (SUSY) models with i?-parity conservation. In the SUSY 
standard model, the lightest neutralino, composed of superpartners of the gauge and Higgs 
bosons, is usually the LSP. From both the theoretical and experimental/observational points 
of view, the wino/higgsino DM is interesting, and it has been studied extensively so far 
. In this paper, we will concentrate on scenarios of the wino/higgsino DM production. 

The present DM abundance is 

n^^Mh'^ = 0.11 ±0.01, (1) 

from the latest WMAP three year data jj], where h is the present Hubble parameter in units 
of lOOkm/s/Mpc. The relic abundance of the LSP must fall in this range to account for the 
DM abundance. The relic LSP abundance depends on the thermal history of the universe. 
Once the universe is reheated up to high temperature after inflation, the LSP is thermally 
produced by particle scatterings and reaches thermal equilibrium. However, it is known that 
the thermal processes cannot provide a sufficient amount of the wino/higgsino LSP with a 
mass, = O(IOO) GeV, because of the large annihilation rates. Therefore, non-thermal 
production must take place below the decoupling temperature, ~ m^/20, to reahze the 
observed DM abundance. 

The non-thermal DM production at such a low temperature may be realized by a late- 
time decay of a scalar field, (p, with a mass m^. Renormalizable interactions generically 
induce a too rapid decay into the standard model (SM) particles and their superpartners 
unless the couplings are suppressed by some symmetries or mechanisms. Thus the decay 
should proceed via non-renormalizable operators suppressed by a large cut-off scale. In this 
paper, we focus on the D = 5 operators. We will show that the branching ratio of the scalar 
decay into the SM superparticles is comparable to that into the SM particles. 

A modulus is one of such late-time decaying fields 0, because the interactions of the 
modulus field are suppressed by the Planck scale, Mp = 2.4 x 10^^ GeV. The decay tem- 
perature of the modulus field is likely lower than the decoupling temperature, leading to 
the non-thermal production of the wino/higgsino LSP. However, according to the recent 
works 0, 3, y, y, Q , the late-time scalar decay is generically plagued with the gravitino 
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overproduction. Thus the cosmological scenario with a modulus field is strongly constrained, 
for example, from the big-bang nucleosynthesis (BBN). 

We point out that a natural source of the non-thermal DM production is related to 
new physics at an intermediate scale, M <^ Mp. The gravitino overproduction problem is 
then relaxed. This is because the total decay rate is enhanced as F ~ m^/M^, reducing 
the branching ratio of the gravitino production. Note also that, since a light scalar field 
becomes acceptable, the decay into the gravitinos may be kinematically forbidden. In the 
following sections we will show that a right amount of the wino/higgsino DM is obtained by 
the scalar decay at an intermediate scale. We propose the saxion field, which is the bosonic 
SUSY partner of the axion, as a natural candidate for such a scalar field. 

Before closing the introduction, let us clarify differences from the works in the past, 
which discussed the non-thermal DM production from late-time decaying scalar fields. The 
decay from a heavy scalar with Planck-suppressed couplings, i.e., a modulus field, has been 
ex„ in Ref. ^. A.O. Ret. Q ana,3.ed a .n,U. .bjec. w.h a ,owe. cu. 

off, although they regarded the branching ratios of the superparticle as a free parameter. In 
this paper, we explicitly show that the branching ratio of the SM superparticle production 
is generically large, and that the effective number of the LSP produced by one scalar decay, 
A^;^, is close to unity. This has a great impact on the non-thermal DM production scenarios. 

The rest of the paper is organized as follows. In Sec.|n]we estimate the branching fraction 
of the SM (super)particle production for the scalar decay through the non-renormalizable 
operators of D = 5. We reexamine the conventional scenario using the modulus decay for 
the non-thermal DM production in Sec. IIIII It is shown that the saxion decay can account 
for the present DM abundance in Sec. IIVI The last section is devoted to conclusion. 



II. HEAVY SCALAR DECAY VIA D = 5 OPERATORS 

In this section we consider a heavy scalar decay into the SM particles and their superpart- 
ners, assuming the scalar mass m,^ much larger than the weak scale. Since those particles 
with R-parity odd eventually produce at least one LSP per decay, the production rate of the 
superparticles is crucial to determine the relic LSP abundance. In the following, we assume 
that the scalar field is R-parity even and singlet under the SM gauge groups. The purpose 
of this section is to estimate the decay rates of the scalar field into the SM (super)particles 
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via the non-renormalizable operators of D = 5. Although such operators have been studied 
in Ref. j^, they overlooked some important interactions of the superparticle production. 
Below we investigate the operators one by one. 

(I) Let us first consider the scalar decay into the gauge boson and gaugino. Such an 
interaction is effectively given by the dilatonic coupling with a cut-off scale M, 

C = ^ j dH (t>W^''^W^''\ (2) 

where Xq is a numerical coefficient and W^'^^ is the supersymmetric field strength of the SM 
gauge supermultiplet. Including the supergravity effects, we obtain 
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+Lo/^{[g%^ + G%cP^) X^^^VrX^^^ + {g%cP^ + X^^^VlX^^^} 



4" l\- <pr • - J ■ ■ ■ Jt- • ■ \- 0T ■ - (pr J ■ ■ ■ ^- ■ jj J (3) 

where the real and imaginary components of the scalar field are defined as = {(pR+i(pi) / V^. 
Although there are the other contributions from the gaugino kinetic terms, the resultant 
decay rate through them is proportional to the gaugino mass squared, m\, after applying 
the equation of motion. Therefore they are neglected unless is close to (but larger than) 
2mx- The subindex, i, attached to the total Kahler potential, G = K + \n iVTp, represents 
a derivative with respect to the field, i, while the superscript is obtained by = g^^'Gj*, 
where g^^* is the inverse of the Kahler metric gij*. The coefficient e*^/^ is equal to the 
gravitino mass at the vacuum, 7/13/2 = (e*^/^). 

The coefficients of the interactions with the gauginos are 

= (g'^^G^)^ = g^^G-^-^ + {g^%G-^, (4) 
G\ = {g^'^G^)^ = 1 + {g'^%G^, (5) 

where we have assumed g^j* = for i 7^ j, for simplicity. If 7/13/2, we can see G^ ^ + 
G'^^ = 0(1) barring cancellations. In fact, the minimal Kahler potential without the SUSY 
mass term (i.e., G^^i, = 0) gives G'^ ^+G'^^ = 1. On the other hand, if the mass is enhanced by 
the SUSY mass, namely ~ e^^'^\G^^\ ^ "^3/25 the coefficient is enhanced as G'^ ^ + G'^^ = 
O^rrifj)/ 1713/2). However if the scalar mass m<^ is enhanced due to the non-SUSY effects, e.g., 
6K = \(f)\^\z\yM^ with the SUSY breaking field z, the coupling is given by G'^^+G"^^ = 0(1). 
Then the decay rates into the gauge multiplets are evaluated as 



I ^\ I ^ TTZ 

Tg = T{(j) gauge boson) ~ T{(f) ^ gaugino) ~ ^g^—lT^i (6) 



either if ~ or if 3> m3/2 due to the SUSY mass term. Here Ng is the number of 
the possible final states, and Ng = 12 for the SM gauge groups, SU{3)c x SU{2)l x U{1)y- 
We find that the production rate of the gaugino is comparable to that of the gauge boson, and 
is not suppressed by the gaugino mass. On the other hand, if the scalar mass is smaller than 
m3/2, the gaugino production rate is estimated as F ~ m^i^m^/ M"^ up to a numerical factor, 
while the gauge boson production rate is the same as Similarly, if is enhanced due 
to the non-SUSY mass term as ^ ^123/2, the gaugino production rate is F ~ m\i2m^/ M'^ , 
which is suppressed compared to the gauge boson production rate F ~ m^/M^. 

(II) The scalar field can couple with the matter fields in the Kahler potential. When the 
Kahler potential includes such a non-renormalizable term as, K = (pQQ'yM + h.c, in which 
the matter fields Q and Q' appear with the opposite chiralities, we can show that the decay 
rate into the SM (super)particle is proportional to powers of the masses of the matter fields. 

Let us first consider the matter scalar production rate. The first contribution comes 
from the kinetic term of the scalar fields: Ck = gij*dn(j)^d'^(f)^* . Then the decay rate is 
proportional to the fourth power of the mass of the matter scalar after using the equation of 
motion, since the interaction term becomes Ck = —{4'/M) Q{d^Q') up to a total derivative 
term, where Q is the scalar component of the matter chiral multiplet Q [5J . Thus the decay 
rate is estimated as Tk ~ {niQi /m^Y x m^/M^. 

The other couplings arise from the scalar potential. In particular, the decay rate may 
be enhanced due to the supergravity effects when m3/2 is large. Since non-renormalizable 
operators violate the conformal symmetry, the operator receives corrections of order at 
the leading level. Thus we obtain Cy ~ J!''^l/2^QQ"^ + though the coefficient depends 
on details of the models. The decay rate then becomes Ty ~ (^3/2 /""^(/i)^ x m^/M'^. 

The above supergravity effects seem to give a correction to the matter scalar mass, Suiq ~ 
"^3/2' tiy taking a vacuum expectation value (VEV), (0) ~ M. It should be noted, however, 
that this estimate is too naive and actually overestimates the correction when m3/2 is larger 
than the weak scale. This is because there exists a cancellation among the contributions 
to the soft scalar mass if \G^\ <^ (0). This cancellation can be understood in terms of the 
Kahler invariance. To be concrete, let us consider the minimal Kahler potential with the 
coupling 6K = (f)QQ"^ /M + h.c. Since the above contribution in 0(171^2) arises from K^W 
of the scalar potential, it is always accompanied by Wtf,. Thus when the F-term of (f) is 
suppressed, the correction of ©(mg/g) will cancel with that from after taking the VEV. 
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Lastly, the production rates of the matter fermions, ipQ^ ^ire chirally suppressed. The 
decay rate is estimated as T^^ ~ [m^^^ / m^Y x m^/M^, which is proportional to the fermion 
mass squared. 

(Ill) The heavy scalar field can decay into the Higgs supermultiplets via the interaction: 

^ = J^l d^O(i)^H^H, + h.c., (7) 

though such a coupling is related to the /x-problem and it is quite mo del- dependent whether 
or not such a coupling exists. We notice that the Higgs chiral supermultiplets appear with 
the same chirality. Expanding the interaction ((Tj), we obtain 
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id^<P^)H^H, + e^/^GV^^fiif J + h.c. (8) 



In addition, the scalar trilinear coupling (^^HuHd arises from the scalar potential; the decay 
rate through the supergravity corrections is given by ~ |A//pm3^2"^0/^^- We neglect the 
interaction terms whose coupling is suppressed by the higgsino mass. In particular, they 
include the U{1) connection term in the covariant derivative of the higgsino kinetic term, 
and that obtained by expanding e*^/^ in the higgsino mass term. These are generically 
small in the mass eigenstate basis by considering the interference with the SUSY breaking 
field BQ. 

Thus, either if ~ m3/2 or if m,^ ^ Tnz/2 due to the SUSY mass term, the decay rates 
are given by 

I \ 12 3 

r(</>^M) ^ V{^^HH) ^ (9) 

for the Higgs multiplets which are doublet under the SU{2)l symmetry. We find that the 
decay rates are comparable to Vg if Ac ~ A//, and that the higgsino production rate is the 
same order of the Higgs production rate. In contrast, when the scalar mass is enhanced by 
the non-SUSY mass term as ^ "^3/2, the higgsino production is suppressed compared 
to the decay into the Higgs bosons by O {171^2 /''tT''^)- 

(IV) The following interactions in the superpotential permit the three-body decay. 
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j d^e(t)HJ'^Q + h.c., (10) 



where Hu,T,Q are the chiral supermultiplet of the up-type Higgs, right-handed top quark 
and left-handed quark doublet of the third generation, respectively. Then the interactions 



become 



<pHj.q + <i)t*HuRq + <i)tHuLq +h.c., (11) 



in the global SUSY Lagrangian. Although the decay rates do not receive chirality suppres- 
sions, the decay is a three-body process and the phase space is small jsf: 

I \ 12 3 

T{ct>-^Htq) ^ r(0->iJtg) ^ 1(0 ^ i^tg) y^^, (12) 

for ^ w^soft• Note that couplings with the lighter quarks and leptons must be suppressed, 
since the corresponding Yukawa coupling constants need fine-tunings by taking account of 
the scalar VEV, (0) ~ M. 

Additionally, the scalar couplings receive supergravity corrections. The decay rate 
through the correction which is of order m^j2 for the operator (j)HuU*q is suppressed both by 
the phase space and by the gravitino mass squared. If < m3/2, the decay rate becomes 
larger than ()12j) . At first sight, the scalar VEV seems to induce a correction of an order 7713/2 
to the trilinear scalar coupling, HuU*q. However, it is actually suppressed if \G^\ <^ (0), 
since the scalar contributes to the A-parameter through the auxiliary component due 
to the Kahler invariance. 

From the analyses of (I) - (IV), either if ~ '"73/2 or if ^ 7773/2 due to the SUSY 
mass ^ 1, almost the same amount of the SM superparticles is produced as that of the 
SM particles: 

Br(0 — i> SM particles) ~ Br(0 SM superparticles), (13) 
with the total decay rate 
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r m, , 

(14) 
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where c is a numerical coefficient, and given by c ~ Ng\^ + \^fj, for the dilatonic coupling and 
the interaction with the Higgs supermultiplets. The other decay processes are suppressed 
either by the soft masses or by the phase space. On the other hand, when the scalar 
mass is enhanced by the non-SUSY mass term as 777^ ^ ^3/2? the above superparticle 
production rate becomes suppressed compared to that of the SM particle. Note that the 
single-gravitino production induce a serious problem in this will be discussed in the 



next section. If 7/13/2 > rn^, the supergravity corrections significantly affect the decay rates, 
and the branching ratio of the superparticle production depends on the non-renormahzable 
operators. We emphasize again that the branching ratio of the superparticle production 
is generically sizable, Br{(f) —>■ SM superparticles) = 0(0. 1), either when ~ ""^3/2 or if 
^ m2,/2 due to the SUSY mass term. 



III. DIFFICULTIES IN MODULUS DECAY 



It has been known that the wino/higgsino LSP with a mass of O(IOO) GeV cannot be 
produced sufficiently in the thermal bath because of its large annihilation cross section. 
We need therefore the non-thermal process below the freeze-out temperature in order to 
explain the observed DM abundance. In this section let us reexamine the non-thermal DM 
production from the modulus decay. 

The modulus field corresponds to a dynamical degree of freedom associated with a fiat 
direction of the scalar potential. Since the SUSY breaking effects stabilize the potential, 
the modulus field naturally obtains a mass around the gravitino mass, 1713/2- In addition, 
the modulus mass can be even larger with sizable non-perturbative corrections. The mod- 
ulus decay is generally induced by the non-renormalizable operators, whose couplings are 
suppressed by the Planck scale, namely M = Mp in the previous section. Such a light and 
late-decaying scalar field is likely to dominate the energy of the universe before its decay. 
Then the subsequent thermal history depends on the decay processes of the modulus field, 
that is, the decay temperature and the branching fractions of the decay products. The 
modulus mass is required to be ^ 100 TeV ^ from the cosmological BBN bound. 

If the decay takes place after the freeze-out of the LSP, the relic LSP abundance is esti- 
mated in terms of the following parameters: the annihilation cross section of the LSP, {cv), 
the modulus decay temperature, T^, and the branching fraction of the superparticle produc- 
tion, which provides an effective number of the LSP from one modulus decay, N^^. According 
to the previous section, the branching fraction of the LSP production is generically quite 
large. Then the resultant LSP abundance after the pair annihilation becomes independent 



One might consider the modulus decay via the D — 6 operators. However it is unhkely that the decay 
takes place before the BBN starts. 
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of the initial abundance, and given by 
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^ 45 1 

V Sn^g^d {(^v) MpTd 

or equivalently, the rehc LSP density is 

where is the LSP mass. In the case of the wino LSP, the annihilation cross section is 
estimated to be e.g. {av) ~ 3.3 x 10^'' GeV^^ for the wino mass M2 = 100 GeV, and that 
of the higgsino LSP is roughly the same order of the magnitude (see Eqs. (jH^ and ((221) )■ 
Since the annihilation cross section decreases for the larger LSP mass, it is likely smaller 
than 0(10^^) GeV~^. Thus the decay temperature must satisfy ^ 100 MeV in order to 
obtain a right amount of the relic DM abundance (^. 

Using the total decay rate (jl4|l with M = Mp, the decay temperature is given by 



3/2 



where c = 0(1). To satisfy ^ 100 MeV, the modulus mass should be large enough: 

> 10^-^ GeV. (18) 

The bound on the modulus mass becomes severer for smaller c and the heavier wino/higgsino. 
It should be stressed again that one of the crucial points to derive the above constraint is an 
effective number of the LSP from the modulus decay. As long as the decay proceeds through 
the non-renormalizable operators of D = 5, the branching ratio of the LSP production is 
generically large, and close to unity. This is another difficulty of the modulus decay, and let 
us call it as "the moduli-induced LSP problem" . 

The mass of the modulus field depends on the mechanism of the potential stabilization. 
A modulus field as heavy as ^ 10®^^ GeV may be realized in the scenarios of the mixed 
modulus-anomaly/KKLT mediation and in the racetrack setups 13[. However, such a 
heavy modulus generically suffers from the gravitino overproduction [6|. Since the conformal 



anomaly mediates the SUSY 
related to the gravitino mass 



Dreaking effects to the visible sector and the correction is 
ll, the gravitino mass is favored to be ^ 100 TeV, leading 



to the relation > 2m^j2- The production rate of the goldstino, which is a longitudinal 



9 



'l-^ 

r(<l> 2V'3/2) ^ — (19) 



c„mpo„e„t of the g^avitino. can be enhanced in the mode , 2^.,_ in high energy Hmit 

^ ^3/2- Actually, the decay rate is obtained as p 

T(eff)|2 ^5 

^3/2 ^^P 

where $ denotes the mass eigenstate (mainly composed of 0), and Q'^^^^ is the effective 
auxiliary field of $. The auxiliary component of the scalar field crucially depends on the 
structure of the SUSY breaking sector as was pointed out in Ref. jol- Since both the modulus 
and SUSY breaking fields are not generally protected by any symmetries at the vacuum, 
they couple with each other, and the rate becomes sizable In fact, the effective auxiliary 
field can be as large as Q^^^ = K3/2"^3/2/'m(^ with K3/2 = 0(1) 10]. Thus the branching 
ratio of the gravitino production is -B3/2 = 0(0.01 — 0.1). Such models with a large gravitino 
production is severely constrained by the cosmology, which we call "the moduli-induced 
gravitino problem" ^. 

As a final remark of this section, let us mention the case of a modulus mass coming 
from the non-SUSY mass which satisfies ^ ^3/2- In this case, a mass of the fermionic 
component of , is hUely much smaller than ra„ so the single-gravitino product.on starts to 
dominate over the pair production. The decay rate is jl3] 

r(^-W'3/.)-j^^ (20) 

for 7710 ^ fn-ip^, "^3/2, which dominates over the other decay processes. It is noticed that the 
rate is independent of the VEV of the effective auxiliary field Q^^^^^ . Thus such modulus field 
generally suffers from the production of too much gravitinos. 

In summary, the modulus decay is disfavored as a source of the non-thermal DM pro- 
duction because of the above "the moduli-induced gravitino/LSP problems." In the next 
section, we propose a natural solution which is free from these problems and also explains 
the relic wino/higgsino DM abundance. 



^ Even if the modulus does not couple to the SUSY breaking sector, the branching fraction of the gravitino 
production is still large. If the scalar mass of the SUSY breaking field, m^, is smaller than as 
B3/2 = 0(10^"' — 10~'^), while -B3/2 = 0(0.01 — 0.1) if > m^. Such large branching ratio causes serious 
cosmological problem [l^ . 
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IV. DARK MATTER FROM SAXION DECAY 



The difficulties of the modulus decay in producing a right amount of DM are associated 
with the fact that the interactions are suppressed by the Planck scale, as shown in the 
previous section. The modulus mass cannot be smaller than 100 TeV to decay before BBN, 
which makes the modulus decay into the gravitinos almost inevitable. One of the attractive 
ways to get around the problem is to postulate a cut-off scale smaller than the Planck scale; 
the branching ratio of the gravitino production can be then suppressed, or may even be 
kinematically forbidden. Thus an intermediate-scale physics is likely to play an important 
role in the non-thermal wino/higgsino DM production. 
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An intermediate-scale physics naturally arises in the Pecci-Quinn (PQ) mechanism 
IitI . Iisl . Q]. The strong CP problem is one of the profound puzzles in SM, and the PQ 
mechanism using a pseudo Nambu-Goldstone boson called axion provides an elegant so- 
lution to the problem. The axion field has a fiat potential at a perturbative level, and 
acquires a finite mass from non-perturbative QCD effects via the anomaly. The axion rolls 
down to the potential minimum at which the strong CP phase vanishes. From the cos- 
mological/astrophysical observations, the axion decay constant, F4, is required to satisfy. 



> 



10^ GeV from SN1987a 
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2l|, and Fa < lO^^-^^ GeV by the axion-overclosure limit ^, 



3, 



24|. 



although the latter depends on the cosmological scenario |22l . 

In the SUSY theories, the axion field is extended to the axion supermultiplet, A. The 
multiplet contains two bosonic real fields, which are axion and saxion, and one chiral fermion 
called axino. The scalar component of the axion multiplet is represented as 

A = -^(a + za), (21) 

where a and a are saxion and axion fields, respectively. The saxion behaves like a modulus 
in the evolution of the universe. Since the saxion potential as well is fiat at a perturbative 
level in the SUSY limit, the field is stabilized by the SUSY breaking (and perhaps non- 
perturbative corrections), and it acquires a mass around the gravitino mass, rricr ~ ?tJ'3/2 
The initial position of the saxion generically deviates from the potential minimum in the 



As is well known, the axion itself is a DM candidate. However the contribution depends on the initial 
condition of the axion field and the subsequent thermal history. In the following, we will focus on the 
abundance of the wino/higgsino LSP. 

Precisely speaking, the saxion mass is model-dependent. For instance, in the heavy gravitino setups, if 

11 



early universe, and it starts to oscillate when the Hubble parameter becomes comparable to 
the saxion mass. Then the energy fraction of the saxion can be sizable, and large number 
of the superparticles may be produced at the decay 

It is a strength of the couplings with the SM sector that significantly differs from the 
modulus field. The axion supermultiplet interacts with the SM (super)particles with the 
couplings suppressed by the axion decay constant, F4. Since -F4 must be much smaller than 
the Planck scale because of the cosmological requirement, the saxion decay rate is much 
larger than that of the modulus field with the same mass. In particular, the axion multiplet 
A necessarily has a dilatonic coupling: 

C = - [ d^6^^Aw^-')W^''\ (22) 

where a is a gauge coupling constant, and F4 the axion decay constant. Here we adopt the 
hadronic axion model h8| for simplicity, and a brief comment on the case of the DFSZ axion 
model [12(1 will be given later. The dilatonic coupling arises from the chiral anomaly of the 
PQ symmetry which is extended to the complex U{1)pq symmetry in the SUSY framework. 
The interaction (j^^ actually reproduces the axion coupling with the SU{3)c gauge bosons, 
providing a solution to the strong CP problem. It also provides the saxion couplings with 
the gauge bosons and gauginos: 



^ ^^p{a)p{a)^,y _ G/2 

871 StcF/ 



(G^a + G^^) A^'^^PlA^'^) + h.c. 



(23) 



where F^"^ and A*^"^ denote the field strength of the gauge boson and the gaugino field, 
respectively. The decay rates of the saxion into the gauge bosons g and gauginos A are given 
byf 

n" - - ^jf. n. ^ AA) . ^\.?^ (24) 



the axion multiplet is sequestered from the SUSY breaking sector, the saxion mass can be smaller than 
the gravitino mass by a loop factor. 

The LSP may be the fermionic superpartner of the axion, axino. However the axino mass depends on the 
models, and we assume that the wino/higgsino is the LSP in the following. 
^ The saxion may have a coupling with the axion, which enables the saxion to decay into two axions. If 
the saxion dominantly decays into the axions, it can be cosmologically problematic |2^. However, it 
is model-dependent whether this process dominates the saxion decay. In fact, such a coupling can be 
suppressed. In the following, we simply assume that the decay mode is suppressed. 
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with 

«-(G^ + G-,)^, (25) 

1 1 I'd 

where we have neglected the sub-leading terms suppressed by {mx/m^Y (m^: the gaugino 
mass). Although we have assumed that the saxion universally couples to the SU{3)c x 
SU{2)l X U{1)y gauge multiplets in estimating the decay rates, this assumption is not 
crucial for the following arguments. This is because the saxion decay into the gauge bosons 
is dominated by that into gluons, to which the QCD axion multiplet necessarily couples. 
Thus the branching ratio of the gaugino production is Br((T AA) ~ + which 

depends on k. As discussed in SecHH G^^ + G^^ is generically 0(1) for < ""^3/25 leading 
to K = 0(1). In fact, for the minimal Kahler potential we have (O"^^ + = 1. Here it 

should be note that the axion multiplet is not allowed to have nonzero SUSY mass, since the 
potential minimum of the axion would be shifted otherwise. On the other hand, if a mass 
hierarchy, mo- ^ ^^3/25 is realized, k can be suppressed as 0(m3/2/mo-). Thus the branching 
ratio of the gaugino production is sizable in the saxion decay, 

Br(a ^ AA) = = O(O.l), (26) 

unless the saxion mass is much larger than m3/2 due to the non-SUSY mass. It should be 
stressed that the branching ratio of the gaugino production is not suppressed by the gaugino 
mass as noted in the previous section, and that the decay rates are enhanced by 0{Mp/ Fj^) 
compared to the modulus case. 

The relic wino/higgsino LSP abundance is determined by solving the Boltzmann equa- 
tions. The wino/higgsino LSP can be produced from thermal scatterings, but it is known 
that the thermal relic abundance is too small to account for the present DM density. As 
will be shown below, the saxion decay can non-thermally produce a large enough number of 
LSPs, due to both large branching ratio of the gaugino production and relatively low decay 
temperature. 

The LSP number density evolves as 

n^ + 3Hn^ = -{av)n\, (27) 

where we have neglected the production of the LSPs from thermal scattering processes by 
assuming that the saxion decay temperature is lower than the decoupling temperature of 
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the LSP, Ty ~ Here (av) is the thermally averaged annihilation cross section of the 

LSP. Let us rewrite Eq. ^T7\i in terms of the LSP abundance, = n^/s, 



dT 



(28) 



where g^^i is an effective number of massless degrees of freedom at the saxion decay, and 
we have approximated that g^d is almost constant during the evolution. By integrating this 



equation, we obtain the simple analytic formula of the relic abundance for T < Td 



3 



YATd)-' + {cyv) Mp{Td - T) 



45 



(29) 



If the initial LSP abundance y^(Trf) satisfies 
YxiTd) » ^ 



45 



^.4 X lo-^^r^ 



{av) MpTd 



■l0.75\^ /lO"^ GeV 



-2^ 



'A 



1 TeV\ 2 

~) VlOi^GeVy 



(30) 



g*d J V (^^) y 

the second term in the parenthesis of Eq. (j29|) dominates over the first term. Then the final 
abundance is approximately given by 

' 1 



45 



(31) 



\87t^g.d {(rv)MpTd 

Therefore the relic abundance is determined only by the annihilation cross section of the 
LSP and the decay temperature of the parent particle, independent of the initial amount of 
the LSP. 

The annihilation cross section, (aw), for the wino LSP is evaluated as 



2TT ^ (2 - xwr 



(32) 



where xw = ^w/^x denotes the mass squared ratio of the W boson and the LSP (wino), 
and g2 is the gauge coupling constant of SU{2)l. For = O(lOf) GeV, (av) is roughly 
estimated as 0(10^'') GeV~^. For the higgsino LSP, it is given by [2a| 



{av) 



9l -2(1 



3271 ^ (2-xvi.)2 ' 
which is smaller than the wino cross section by one order of magnitude. 



(33) 
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Assuming that the saxion dominantly decays into the gauge bosons and gauginos with 
the rates given by Eq. (|^. the decay temperature is given by 

'10.75\^ f \l /lO^^GeVX 



,,.,,xlO'^MeV,ra ^ . (34) 



where ^ is defined as ^ = y(l + |kP)/2. Using Eqs. and (jMj) . the LSP rehc density is 
given by 



10.75\^ fW-^ GeV"^\ / \ /I TeV\' 



as long as the saxion decay provides a sufficient amount of the LSP. 

The initial LSP abundance produced from the saxion decay depends on the saxion den- 
sity at the decay as well as the branching ratio of the gaugino production. Since the saxion 
potential is generally fiat, it can develop a large expectation value during infiation like a 
modulus field. When the Hubble parameter becomes comparable to the saxion mass, it 
starts to oscillate around the potential minimum with an amplitude at. The cosmological 
abundance of the saxion thus depends on the initial displacement from the potential mini- 
mum, which is expected to fall in the range from F4 to Mp. The saxion dominates the the 
energy density of the universe if the following condition is satisfied: 

1 

O-i n: (T, 

> 



V-e{^)\ (36, 



Mp MR. 

where Tp is the reheating temperature, and we have assumed that the reheating is not 
complete when the saxion starts to oscillate. If the saxion dominates the universe, the 
saxion-to-entropy ratio at the decay is given by 

independent of the initial amplitude o"j. The initial LSP abundance is related to as 
Y^(Td) = N^Ya, where is the averaged number of the LSP produced by the decay of one 
saxion field. Since Ny. is close to unity from Eq. ()2fi|l . Y^{Td) easily exceeds the critical value 
Yj^\ and therefore the relic LSP abundance saturates to the value given by (jHT|) . 

Even when the saxion field does not dominate the universe before its decay, the LSP com- 
ing from the saxion decay tends to exceed the critical abundance Y.^'^K Since the reheating 



15 



temperature is then bounded above due to the gravitino problem |27| ^, the saxion is hkely 
to start to oscillate before the reheating process completes. This is the case if Tr < T^-, 
where T^. is defined by 

Here g^:osc counts the effective number of massless degrees of freedom a.t H = m„. The 
saxion abundance is 

°-"f^)(lo^)(f)^(lo4vV- (-) 



Thus even if the saxion does not dominate the universe, (jHn|) is satisfied for ctj ^ F4. So, in 
the following, we adopt Eq. to estimate the relic LSP abundance, assuming that (IHUj) is 
satisfied. Note that it is the large branching ratio into the gauginos that enables the saxion 
to non-thermally produce a right amount of the LSP DM, even if it does not dominate the 
universe. 

In Fig.Q we show the contours of the LSP abundance in (mg-, -F4) plane. We have chosen 
the wino LSP with = 300 GeV and ^ = 1 as representative values. Since the gaugino 
branching ratio is generically 0(1) as long as mn-_> 2mv, the relic LSP abundance is given 



by Eq. (jHSj) . From the observation of SN1987A |20|,|21|, F4 is bounded from below, F4 > 10 



^9 



GeV, while an upper bound, F 



< 10 GeV, is set by taking account of an overclosure limit 



of the relic axion abundance j22,|23[. We find that the observed DM abundance, Q^h"^ ~ 0.1, 
is realized for rrif^ around O(10^~'^) GeV and F4 = 10^°^^^ GeV (below the dashed line). 
Note that the axion overclosure bound can be relaxed by the entropy production associated 
with the saxion decay. In fact, in deriving the upper bound F4 < 10^^ GeV, it is assumed 
that the initial displacement of the axion from the true minimum is 0(1), and that there is 
no entropy production [28j after the axion starts to oscillate at T = Aqcd — 200 MeV. If 
the saxion dominates the universe and decays after the QCD phase transition (of course, but 



before the BBN begins j29|, |30|), the upper bound can be relaxed to ~ 10^^ GeV j2^. Then 
although it depends on the initial VEVs of the axion and saxion fields, the region above the 
dashed line in Fig. ^ can be viable, and the right amount of the wino DM is realized even 



8 Throughout this paper we assume that there is no source for the late-time entropy production [23 other 
than the saxion. 
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FIG. 1: Contours of the wino LSP abundance, ?L^h^ = 0.3, 0.1, 0.03, 0.01, for = 300 GeV. 
The shaded regions are excluded; the LSPs are only thermally produced and therefore too small 
relic density for mg- < 2m^; the upper left triangle region is excluded since the saxion decays after 
the BBN starts, i.e., < 5MeV, while the lower right one corresponds to the decay temperature 
higher than the decoupling temperature of the wino, i.e., > m^/2Q] J^^^^ ^ ^-^^ from the 
LSP overproduction in the light of WMAP; < 10^ GeV from the SN bound. Also the region 
above the dashed line is disfavored by the axion overclosure limit if no entropy production occurs 
after the QCD phase transition. 



for that region ^. 

In the above numerical analyses, we have fixed the wino mass as m-^ = 300 GeV. In fact, 
the larger LSP mass favors smaller for fixed mo-. This can be understood as follows; higher 
saxion decay temperature is needed to realize f^^^^ ~ since larger suppresses the 
annihilation cross section (see Eq. ()32|)). It is interesting to see that for ^ 500 GeV, the 

^ The DM abundance includes the contribution from the relic axion as well as the wino/higgsino LSP. When 
we take into account the axion contribution, the smaller LSP abundance tends to be favored, though the 
result depends on both the initial displacement of the axion and the thermal history of the universe. 
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parameter region to reproduce Q^h'^ ~ 0.1 (calculated by using Eq. overlaps with that 
where the LSP is thermalized, namely, ^ m^/20. This means that the saxion decay takes 
place before the decoupling of the LSP. Then we need to numerically solve the Boltzmann 
equation around the LSP decoupling, taking account of the LSP production from the saxion 



;he relic abundance depends 
we do not go into details 



decay in order to estimate the LSP abundance correctly. Since 
sensitively on the evolution around the decoupling temperature 
in this paper. On the other hand, the lighter LSP may be excluded by the BBN. The ^Li 
can be overproduced due to the energetic hadrons from the DM annihilation. Thus a mass 
of ^ 100 — 200 GeV is disfavored for wino/higgsino DM 8^ . 

So far we have considered the wino LSP. When the LSP is the higgsino, ^^^"^ ~ ^-^ 
requires higher decay temperature. In fact, since the annihilation cross section is lowered by 
~ 1/10, the decay temperature should be raised by about one order of magnitude. Then the 
constant contours of the LSP abundance in Fig. ^ shift to lower F4 by ~ 1/10 for fixed 
Therefore, we find that the observed relic abundance of the higgsino DM can be naturally 
realized by the non-thermal production from the saxion decay, for rria at the weak scale and 

In the above discussions, we have considered the hadronic axion model. Let us now 
comment on the DFSZ axion model. In this case the standard model particles have nonzero 
PQ charges, and the saxion can dominantly decay into the third generation quarks. Then 
the total decay rate of the saxion becomes slightly higher and A^^^ is suppressed by that 
amount. It should be noted, however, that the relic DM abundance is given by Eq. 
as long as is satisfied. Therefore the contours of the LSP abundance in Fig. ^ slightly 
move upward for the DFSZ axion model, but our main conclusion remains virtually intact. 

Lastly let us comment on the gravitino production in the saxion decay. If A has a mixing 
with the SUSY breaking field z, the saxion can decay into a pair of the gravitinos, in addition 
to the single gravitino production. In contrast to the modulus, however, the saxion decay is 
free from the gravitino overproduction problem. This is because the axion decay constant 
is much smaller than the Planck scale. Since the total decay rate is enhanced by lowering 
the cutoff scale, the saxion mass needs not be larger than the gravitino mass to decay 
before BBN, and the gravitino production may be kinematically forbidden. Even when the 
channel is kinematically allowed, the branching ratio of the gravitino production can be 
suppressed. Let us first consider the gravitino pair production rate. Since the coupling of 
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A with the two gravitinos is proportional to -F4 [10], the branching ratio of the gravitino 
production is suppressed by (F^/Mp)'^, namely = 0{Fj^/ MpY . Since F4 is bounded 
above, F4 ^ 10^^ GeV, leading to K3/2 ~ 10^^^ which is small enough to evade the bounds 
from the BBN and LSP overclosure limits [6]. However, the single-gravitino production rate 
actually dominates over the pair production rate, if the saxion mass is larger than m^i2 
due to the non-SUSY mass term. The branching ratio of the single-gravitino production 
is given by 0{{m^/ 171^/2)^ {F_a/ Mp)"^) , and the cosmological bounds are not so severe unless 
^ m3/2. Note that, since the saxion does not necessarily dominate the universe to 
produce a right amount of the wino/higgsino DM, the gravitino abundance from the saxion 
decay may be suppressed. Therefore the gravitino production can be safely neglected in the 
saxion decay. 



V. CONCLUSIONS AND DISCUSSION 



In this paper we have pointed out that an intermediate scale physics is necessary for a 
successful non-thermal production of the wino/higgsino DM, to avoid the moduli-induced 
gravitino/LSP problem recently pointed out in Refs. Q]. This is because gravitinos and/or 
LSPs would be overproduced otherwise. The conventional scenario using the modulus de- 
cay to produce the wino/higgsino DM non-thermally is disfavored from this point of view. 
Instead we have proposed an alternative candidate: the bosonic SUSY partner of the axion, 
i.e., the saxion. It is particularly interesting that the right amount of the DM can be real- 
ized for the saxion mass around the weak scale when we consider the axion decay constant 
within the cosmologically allowed range, F4 = 10^^^^ GeV. Although our analyses focused 
on the case of the saxion, we would like to stress that they are rather generic and can be 
applied to any scalar fields that decay via the non-renormalizable operators of D = 5 with 
an intermediate cut-off scale. 

It is also worth noting that the saxion decay is applicable for a source of the harm- 
less entropy production. As was stressed in this paper, the late-time decay of the scalar 
field generally suffers from the gravitino/LSP overproduction. In contrast, the scalar field 
at an intermediate scale can reheat the universe avoiding overproduction of the grav- 
itino/superparticles. For instance, a scalar field at M <^ Mp may decay before BBN begins, 
while satisfying < 2m3/2,2mLSP to kinematically block the gravitino/LSP production 
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channels. 
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